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abstract Finite iV-unit Langevin models with additive and multiplicative noises have been 
studied with the use of the augmented moment method (AMM) previously proposed by the 
author [H. Hasegawa, Phys. Rev E 67, 041903 (2003)]. Original iV-dimensional stochastic 
equations are transformed to the three-dimensional deterministic equations for means and fluc- 
tuations of local and global variables. Calculated results of our AMM are in good agreement 
with those of direct simulations (DS). We have shown that although the effective strength of the 
additive noise of the iV-unit system is scaled as /3(N) = (3(l)/v~N, it is not the case for multi- 
plicative noise [a(N) ^ oc(V) / 'y/N], where a(N) and f3(N) denote the strength of multiplicative 
and additive noises, respectively, for the size-iV system. It has been pointed out that the naive 
assumption of a(N) = a(l)/\/~N leads to result which violates the central-limit theorem and 
which does not agree with those of DS and AMM. 
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It has been recognized that stochastic (Langevin) equations subjected to additive 
and/or multiplicative noises may be good models for discussing the property of many sys- 
tems not only in physics but also in biology, chemistry, economy and networks. In recent 
years, much attention has been paid to multiplicative noise whose interesting aspects have 
been intensively investigated (for a recent review, see ref. 1: related references therein). 
It has been shown that multiplicative noise induces the phase transition where multi- 
plicative noise creates an ordered state while additive noise generally works to destroy the 
ordered state. The effect of multiplicative noise in Langevin equation has been discussed 
in connection with the non-extensive thermodynamics [21 non-Gaussian distribution 
is realized for Gaussian multiplicative noise jUEJIB]. 

The property of multiplicative noise is different from that of additive noise in some re- 
spects. In order to make our discussion concrete, let's consider the simple, linear Langevin 
model given by 

^fi) = _ A x .(t) + a xi{t) vM + (i = l-N) (1) 

where A denotes the relaxation rate, a and (3 the strengths of multiplicative and additive 
Gaussian white noises, respectively, given by zero-mean random variables T]i (t) and 
with < rji(t) Tjj{t') >=< £i(t) £j(t') >= 5ij5(t — t'). The system given by eq. (1) has been 
studied by using the Fokker-Plank equation (FPE) mainly for N = 1 or iV = oo: for the 
latter case, the mean- field and diffusion approximations are employed PjJ [7j . In order to 
discuss the dynamics of the finite- N Langevin model, we define the global variable given 
by X(t) = (1/N) J2i x i(t), an d sum eq. (1) over i to get 

~~dT~ = X(t) + a m { EMt) ) + n^ m ' { ) 

It is reasonable to replace the additive noise (third) term in eq. (2) by (3(N)£ where £ 
denotes the white noise with 

(3) 

As for the multiplicative noise (second) term in eq. (2), Munoz, Colaiori and Castellano 
claimed in a recent paper jH] that it may be replaced by a(N)Xr) with the white noise r) 
in the weak-noise case with 

^MHlrfH- <4) 

if Xi and rji are uncorrelated. Equations (3) and (4) lead to the Langevin equation for X 
given by 

<1X[t) -Ai(t) + 4=x(i)#) + 4ew. (5) 



dt VN VN 

from which the FPE for the finite- N system is easily obtained. The assumption given by 
eq. (4) naively seems reasonable. It may be, however, not justified because Xj and rji are 
correlated: Xj may be increased (decreased) for Xi > (xj < 0) by a positive rji in eq. (1). 
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The purpose of this letter is to examine the validity of the approximation given by 
eq. (4). We will employ the augmented moment method (AMM) (or dynamical mean- 
field approximation) developed by the author, which has been successfully applied to a 
study of coupled stochastic systems 0. In this approach, we have discussed dynamics 
of coupled stochastic systems in terms of a fairly small number of quantities relevant to 
local and global variables, by transforming original stochastic equations to deterministic 
equations. It is essentially the second-order moment method both for local and global 
variables. In ref. 9, we obtained equations of motions for means, variances and covari- 
ances, by expanding variables around their mean values. We have adopted, in this letter, 
an alternative FPE method to obtain equations of motions of the quantities, in order to 
avoid the difficulty due to the Ito versus Stratonovich calculus inherent for multiplicative 
noise. Numerical results calculated by using our AMM approach are in good agreement 
with direct simulations (DS) for the Langevin equation given by eq. (1). On the contrary, 
the result obtained from the FPE using the approximation given by eq. (4) disagrees with 
those of DS and AMM. 

We have adopted the finite A^-unit general Langevin model given by 

^ = F( Xi ) + a G{x t ) Vi(t) + P £,(*) + lt\t), (i = l-N) (6) 

with 

4 e Ht) = j y EMt)-x i (t)], (7) 

where F(x) and G(x) denote arbitrary functions of x, w the coupling strength, Z = N — l, 
a and (3 denote the strengths of multiplicative and additive noises, respectively, and rji(t) 
and express zero-mean Gaussian white noises with correlations given by 

(Vi(t) Vj(t')) = S^it-t'), (8) 

Ut) &(0> = M(*-*0. (9) 
(»*(*) 6(0) = °- ( 10 ) 

The Fokker-Planck equation for the distribution of p({xi},t) is given by JU] 

+ 1 E ^a{W 2 G{x k f + p({ Xl }, t)}, (11) 

Z k OX k 

where G'(x) = dG(x)/dx, and = 1 and in the Stratonovich and Ito representations, 
respectively. 

When we consider the averaged, global variable X given by 

*(*) = ^2>(*)> (12) 

i 

the Fokker-Planck equation P(X, t) for X is given by 

P(X, t)=J---Ju t dx t p({ Xl }, t)s(x-^j2 x i)- ( 13 ) 
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The k-th moments of local and global variables are defined by 

(x!l) = JlUdxipttxilrfx*, (14) 

(X k ) = JdXP{X)X k . (k — 1,2, ••) (15) 

By using eqs. (11) and (14), we get equations of motions for mean, variance and covariance 
of local variable (xj) given by 

= (^)) + |EK^)-(^)] + ^(G , (x,)G(x J )), (16) 

= ( Xi F( Xj )) + (Xj F( Xi )) + | J2[( x i x k) + (XjXk) - - 

+ ^[(xi&ixjMxj)) + (xjG'(xi)G(xi)}} 

+ [a 2 (G(x 4 ) 2 >+/3 2 ]^, (17) 

Equations of motions of mean and variance of global variable (X) are obtainable by using 
eqs. (12), (16)-(17): 

d{X) 1 \- d(xj) . . 

d ( x2 ) 1 v^y^ d(xj xj) 

- jp\\^t- (19) 

The mean-field approximation employs only eq. (16) jTTj. Equations (16) and (17) are 
adopted for a discussion on the fluctuation- induced phase transition in infinite- N stochas- 
tic systems [T2j. Equations (18) and (19) play a crucial role in discussing finite- N systems, 
as will be shown below. 

In the AMM [9], we define the three quantities, /i, 7 and p, given by 

H = <X>=l^<x,>, (20) 

i 

7 = ^£<(x,-/i) 2 >, (21) 

i 

p = <{X-pf>. (22) 

It is noted that 7 expresses the averaged fluctuations in local variables (xi) while p denotes 
fluctuations in global variable (X). Expanding X{ in eqs. (16)-(19) around the average 
value of \i as 

Xi = fi + 8xi, (23) 

and retaining up to the order of < (5xi) 2 >, we get equations of motions for /x, 7 and p 
given by 

^ = fo + f2l + (^y~ J [go9i + ^(9192 + 9093)1} , (24) 



Tt = 2 /i7+(^)(p-7) + (0+1)(^i 2 + 2^2)« 2 7 + «\ 2 + /3 2 , (25) 
| = 2f lP + ((P + lM + 2g g 2 )a>p + ^f + ^, (26) 

where fa = (l/£\)d e F(p)/dx e andg^ = (1/ £\)d e G(p) / dx e . Original A^- dimensional stochas- 
tic equations given by eq. (6) is transformed to three-dimensional deterministic equations 
given by eqs. (24)- (26). In the limit of a — 0, they reduce to those obtained previously 
for the Langevin model subjected only to additive noises [eqs. (7)-(9) in ref. 9c]. 

It is easy to see that the stationary solutions of 7 and p for w = satisfy the central- 
limit theorem: 

P = 1 (27) 

<* 2 9o+P 2 (2R) 

7 " [-2f 1 -(l + ( f ) )(gl + 2g g 2 )a'Y 1 ) 

It is noted that even for N = 00, local fluctuations exist (7 ^ 0), although global fluctu- 
ations vanish (p = 0). 

Now we consider the result when the approximation given by eq. (4) is employed (such 
result is referred to as the APP hereafter). The APP leads to equations of motion for p 
and p given by 

^ = /0 + /27+ (^^j \3o9i + 3(0102 + 9093)1], (29) 

I = 2/ lP+ (l + 0)( 9 ? + 2 9o92 )^ + ^M + £ (30) 

We note that the third term of eq. (29) and the second term of eq. (30) are different from 
their counterparts in eqs. (24) and (26). The ratio between p and 7 in the stationary 
state becomes 



P APP 1 / -2f 1 -(l + <f>)(g 2 + 2g g 2 )a 2 



1 APP ~ N \—2fi — (1 + 0)(g 2 + 2g Q g 2 )a 2 /N ^ 
1 f-2f 1 -(l + ( f>)(g 2 + 2g g 2 )a 2 \ 



(31) 

(as N ^00) (32) 



N V -2/i J 

which is in contradiction to the central-limit theorem. 

In order to make numerical calculations, we have adopted the linear Langevin model 
given by eq. (1): 

F(x) = -Xx, G(x) = x. (33) 
In the Stratonovich representation, equations of motion given by eqs. (24)- (26) become 

da a 2 a 

it = (34) 

^ = -2A 7 + ^(p-7) + 2a 2 7 + aV + /5 2 , (35) 
- = -2\p + 2ap+ — + -. (36) 
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The stationary values for w = are given by 

fi = 0, 

7 : 

P z 



P 2 



2(X-a 2 ) 
1_ 
N' 



(37) 
(38) 

(39) 



Equation (38) shows that a 2 < X because 7 > 0. 

The probability p(x) for the linear Langevin model (N = 1) in the stationary state is 
given by |H Ej 

1 r . 
Z 



p(x) 



1 — (1 — q)rx* 



l-q 



(40) 



with 



-.B 



Q 



(g-l)r V 2 1- 1 

2X + a 2 
2f3 2 ' 
2A + 3a 2 
2X + a 2 ' 



(41) 

(42) 
(43) 



where B(a,b) is the beta function. By using p(x) given by eqs. (40)-(43) with eq. (14), 
the first and second moments are given by 



< x > 



< x > 



0. 



1 



(q — l)r 
1 



B 



1 _\_ 

2' q-1 
P 2 



2 1 



r(5 - 3q) 2(A - a 2 

which agree with the result given by eqs. (37) and (38). 
If we employ the APP, eqs. (30) yields 



(44) 
(45) 

(46) 





(3 2 /2N 


P APP = 


X-a 2 /N' 




P 2 




2iVA' 


p APP 


1 / A-« 2 


jAPP 


\X-a 2 /N 




1 /A-a 2 \ 




n[ x r 



fas — ► 00) 



fas — > 00) 



(47) 

(48) 
(49) 

(50) 



The value of p APP in eq. (48) is independent of a for — > 00. Equation (50) is in 
contradiction to the central-limit theorem. 
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In order to examine the validity of our AMM approach, we have made direct simu- 
lations for the iV-unit Langevin model given by eq. (1). The Heun method is employed 
with a time step of 0.0001. Figure 1 shows the N dependence of the stationary value of 
p for the three sets of parameters: (1) a = 0.8, (3 = 1.0, (2) a = 0.5, (3 = 1.0, and (3) 
a = 0.5, P = 0.5 with A = 1 and w = 0. Squares, circles and triangles show results of 
direct simulations (DS) for the sets (1), (2) and (3), respectively. We note that p follows 
the relation: p oc N~ x . In the case of w = under consideration, we get p = and 
j(N) = p(N = 1) independent of N. Solid curves show the result of the AMM, which 
are in good agreement with those of DS. In contrast, the result with the APP shown by 
dotted curves disagrees with the result of DS and AMM. For N > 10, p APP for the sets 
(1) and (2) becomes almost independent of a as eq. (48) shows. 

We have tried to obtain the Fokker-Plank equation P(X,t) for the global variable X. 
We assume that the FPE for P(X,t) is expressed by 



3P _ 3 f\„ lin , (t±3^ G . (x)G[x) 



Of 




p 



G (X) 2 + ^f] P 



(51) 



where = 1 and in the Stratonovich and Ito representations, respectively, and A, B 
and C are parameters whose values will be determined such as to yield the correct first 
and second moments, as shown below. By using P(X,t) in eq. (51), we get the equation 
of motion for p and p given by 

^ = fo + hi + ( (0 + f ^ ) \go9i + + goftb], (52) 

-£ = 2f 1 p+( ( f ) + A + B)(g 2 1 +2g g 2 )a 2 p + Ba 2 g 2 + C(3 2 . (53) 

A comparison between eqs. (52) and (53) with eqs. (24) and (26) leads to 

<J) + A = 0, (54) 
<j) + A + B = + 1, (55) 

B = C = ± (56) 

Unfortunately, conditions (54)- (56) are satisfied only for N — 1 with A = and B = C = 
1, but have no solutions for N ^ 1. This implies that the FPE for the global variable X 
in iV-unit Langevin model is not expressed by the form given by eq. (51) with A oc N~ s 
for iV > 1 (5: an index). It is probable that J2i x i{t) Viif) 111 the second term of eq. (2) 
may yield non-Gaussian multiplicative noise because the distribution of £j(t) does not 
follow Gaussian. 

To summarize, we have studied finite A-unit Langevin model including additive and 
multiplicative noises by using AMM jU], which has been reformulated with the use of 
FPE. It has been pointed out that the scaling assumption given by eq. (4) adopted in 
a recent paper (ref. 4) leads to results violating the central-limit theorem [eq. 32] and 
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disagreeing with those of DS and AMM. The scaling relation of the effective strength 
against N of multiplicative noise is different from that of additive noise. Our AMM may 
be applied to the general Langevin model with arbitrary forms of F(x) and G(x), and 
also to multi-variable stochastic models like FitzHugh-Nagumo neuronal model. 

A disadvantage of our AMM is that its applicability is limited to weak-noise cases. On 
the contrary, an advantage of the AMM is that we can easily discuss dynamical property 
of the finite iV-unit Langevin system by solving the three-dimensional ordinary differential 
equations. In contrast, within direct simulation and the FPE approach we have to solve 
the iV-dimensional stochastic Langevin equations and the (N + l)-dimensional partial 
differential equations, respectively. Although the discussion presented in this letter is 
confined to the stationary solution of equations of motions given by eqs. (24)- (26) [or 
eqs. (34)-(36)], it is possible to discuss the dynamical property of the coupled Langevin 
model by solving them. Actually in our previous papers 0, we have studied the N- 
dependent synchronization in networks described by Langevin, FitzHugh-Nagumo and 
Hodgkin-Huxley models subjected only to additive noises with global, local or small- 
world couplings (with and without transmission delays). It is interesting to make such 
calculations by including also multiplicative noises, which is left our future study. 
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Figure 1: (Color online) The N dependence of the stationary p, fluctuations in the global 
variable of X, for the three sets of parameters: (1) a = 0.8, (3 = 1.0, (2) a = 0.5, (3 = 1.0, 
and (3) a = 0.5, j3 = 0.5 with A = 1 and w — 0, calculated by direct simulations with 
1000 trials (DS: squares, circles and triangles for (1), (2) and (3), respectively), the AMM 
(solid curves), and the APP (dotted curves), bars denoting the root-mean-square values 
of DS results. 
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